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We consider a Ginzburg-Landau model for superconduc- 
tivity with a Chern-Simons term added. The flow diagram 
contains two charged fixed points corresponding to the tri- 
critical and infrared stable fixed points. The topological cou- 
pling controls the fixed point structure and eventually the 
region of first order transitions disappears. We compute the 
critical exponents as a function of the topological coupling. 
We obtain that the value of the v exponent does not vary 
very much from the XY value, uxy = 0.67. This shows that 
the Chern-Simons term does not affect considerably the XY 
scaling of superconductors. We discuss briefly the possible 
phenomenological applications of this model. 
Pacs numbers: 74.20De, ff.fOHi. 



I. INTRODUCTION 

The Ginzburg-Landau model (GL) has been intro- 
duced almost half a century ago m as a phenomenologi- 
cal model for superfluidity and superconductivity. In its 
one component order parameter version it has been used 
with remarkable success as a statistical mechanics model 
for the critical phenomena of systems lying in the same 
universality class as the Ising model |p[ . In its N compo- 
nent version coupled to Abelian gauge fields it has been 
used as a model for superconductivity and liquid crystals 
||. However, in this last situation the e-expansion, which 
works very well in the non-gauged version, seems to be 
insufficient to describe unambiguously the critical prop- 
erties of the model [QPJ. Only in the large N limit the 
e-expansion gives consistent results [Q . The trouble is the 
absence of a second order (infrared stable) fixed point in 
the flow diagram. This result is physically correct only 
in the extreme type I regime u/e 2 << 1 (here u is the 
quartic scalar self-coupling and e is the charge) , where we 
expect a weak first order phase transition. This regime is 
also well described by the fluctuation corrected mean field 
analysis of Halperin et al. Q . The weak first order tran- 
sition has been probed experimentally in certain classes 
of liquid crystals, where essentially the same GL model 
holds [||. For the extreme type II region (u/e 2 » 1), 
a second order fixed point is expected. Indeed, this re- 
sult follows from numerical studies performed in a lattice 
dual GL model Q|. Therefore, the prediction of a first 
order transition even in the type II regime seems to be 
an artifact of the e-expansion. More recent works sup- 
port also this point of view [||,[l0 11 Also, it has 



been shown that the renormalization group (RG) in a 
fixed dimension approach is more appropriate J|,|l0|,[ll]] . 
The flow diagram in the u — e 2 plane exhibits in gen- 
eral four fixed points: Gaussian and supcrfluid (or XY), 
both uncharged, and the tricritical and superconducting 
which are charged fixed points. The Gaussian fixed point 
is trivial and describes the mean field critical behavior of 
a 0(2) model (we shall work with N = 2). The super- 
fluid fixed point describes the A-transition in ffe 4 and 
lies in the same class of universality as the XY model. 
The tricritical fixed point is over a line, called tricritical 
line, which separates the regions of first and second order 
phase transition. This fixed point is attractive along the 
tricritical line but repulsive in the direction nearly paral- 
lel to the u-axis. Finally, the superconducting fixed point 
is a charged infrared stable fixed point and describes a 
second order superconducting phase transition. Berger- 
hoff et al. |l(]] obtained this flow diagram using Wilson's 
RG [ p"2[ in a non-pcrturbative version called exact renor- 
malization group Jl3| . They used the background field 
formalism to control gauge invariance which is in prin- 
ciple violated due to the presence of the cutoff (further 
discussions on this subtle point in Wilson's RG can be 
found in refs. 14 1). This flow diagram has been also 
obtained by Herbut and Tesanovic [ pd[ using a simpler 
method. They performed a 1-loop calculation in a fixed 
dimension approach. 

Recently a flow diagram with qualitatively the same 
structure has been also obtained for a GL model with a 
topological Chern-Simons (CS) term by Malbouisson et 
al. . Their analysis was performed using Wilson's RG 
in perturbative form, which is not manifestly gauge in- 
variant due to the cutoffed integrals. Due to the presence 
of a CS term, an intrinsically three dimensional object, 
their calculations were performed in d — 3 and N = 2, re- 
sulting therefore in an uncontrolled approximation since 
there is no small parameter as e or l/N. The same type 
of model has been considered earlier by Kleinert and 
Schakel using a different scaling. They performed 
a 1-loop calculation of critical exponents. However, their 
scaling did not allows for a consistent zero CS mass limit 
since Feynman graphs were evaluated at zero external 
momenta. The reason comes from the graph shown in 
Fig.l. When evaluated at zero external momenta this 
graph gives zero due to the structure of the CS gauge 
field propagator. On the other hand, in the same scaling 
this graph is infrared divergent if no CS mass is included 
in the action. Thus, it is not legitimate to perform the 
zero CS mass limit in the scaling considered in Ujn as ob- 
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served by the authors themselves. There are many other 
RG studies of bosonic CS models in the literature but the 
F 2 term is almost always absent [H. The presence of 
such a term is crucial in order to obtain the flow diagram 
of ref. |16|. Moreover, it is desirable to recover the usual 
GL model in the limit of zero CS mass. 

In thi s pape r we consider further the topological model 
of refs. pL 6|Jl 7f| , performing the calculation at the critical 
point |2|. We consider two different approximations. In 
a first step (section II) we perform a 1-loop calculation 
of the RG functions assuming that the same scale holds 
for both the order parameter and the gauge field. In this 
context we find the following main features: (1) For the 
CS coupling smaller than a certain critical value there 
are no charged fixed points; (2) There exists a interval 
of CS couplings such that two charged fixed points are 
found, corresponding respectively to the tricritical and 
second order (infrared stable) fixed points. In this in- 
terval it is possible to find respectable values for the v 
exponent but not for the r\ exponent, eventually violating 
the scaling relations; (3) For larger CS coupling, outside 
the interval mentioned in (2), the region of first order 
behavior is lost since the tricritical fixed point assumes 
an unphysical value. The second step (section III) con- 
sists in improving upon the 1-loop result of section II by 
distinguishing the scales of the order parameter and the 
gauge field. In this approximation we follow an idea of 
Herbut and Tesanovic in order to relate both scales and 
obtain in this way the RG flow. In this case we have 
the tricritical and second order fixed points even in the 
limit of zero CS coupling, consistent with the fixed point 
structure of the conventional GL model. However, once 
again the tricritical fixed point becomes unphysical for 
the CS coupling larger than a certain critical value. In 
this case we obtain that all the critical exponents have 
respectable physical values. Moreover, the v exponent as 
a function of the topological coupling does not deviate 
very much from XY scaling. Finally, in section IV we 
discuss our results and the possible applications. 
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where go = Att^lqOq. 

Now we write the renormalized action as a sum of 
S' + 5S where S' is the same as S with the bare quan- 
tities replaced by renormalized ones (in our notation it 
corresponds to drop the zeroes). SS is the counterterm 
part which is given by 
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where the renormalized fields are given by Aq = \/Z~aA 
and ipo = \J Z^ip. We shall perform our calculations at 
the critical point and therefore r = 0. Also, we choose 
ro in such a way as to cancels the tadpole graphs. The 
renormalization conditions for the renormalized vertex 
functions are given by 
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II. MODEL AND RG RESULTS 



Our starting point is the following action: 
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where the subindex denotes bare quantities. The above 
action is a standard GL model with a CS term added. 
Sgf is the gauge fixing term which is given by 



■%.,- = ! '/ :i -'-^;v-.i;,r. 



(2) 



with «o being the bare gauge fixing parameter. The bare 
propagator for the gauge field is given by 



where u is a dimensionless coupling and S.P. denotes the 
symmetrical point which we take as being given by 



k a - k b = (iSab - 



(11) 



P^p (k) and (k) are the transverse and longitudinal 
projections given respectively by Pj?(fc) = 5y — kikj/k 2 
and Pj? (k) — kikj/k 2 . Note that in the above renor- 
malization conditions the gauge coupling and the scalar 
coupling are fixed at the same momentum scale. 

The renormalized couplings will depend on the momen- 
tum scale p and the beta functions are defined through 
derivatives of the couplings with respect to log p. Note 
that in order to preserve the 0(2) symmetry of the four 
point function, we must consider the Feynman graphs 
with all the external momenta incoming at the vertices. 
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If we proceed otherwise, that is, if we choose a convention 
which two external lines are incoming while the other two 
are outcoming, then the corresponding four point func- 
tion is not symmetric when using the above symmetrical 
point. This is easily seen by performing a 1-loop calcula- 
tion in the uncharged model with O(N) symmetry. The 

resulting four point function T^l s (the subindices are 
color indices) is not proportional to the O(N) symmet- 
ric tensor (SapS^s + d al 5ps + S a sS/3j)/3 if the graphs are 
not evaluated with the convention that all the external 
momenta are incoming. 

The Ward identies imply the exact relations Z e Z^ = 1, 
e 2 = CqZa and ao = olZa, where Z e is the charge renor- 
malization. Let us define the following dimensionless 
gauge couplings, e 2 = e 2 /p, 9 — 9/p, f = Anfie 2 and 
g = Att^iO. The flow equations are given up to 1-loop 
order by 
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where t]^ is the anomalous dimension for the scalar field 
defined by 
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while i]a is the anomalous dimension for the gauge field 
which is defined by 



Va =P- 
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u> is a function of g which will be written later, together 
with the explicit expressions of the corresponding anoma- 
lous dimensions. Before doing this, let us discuss some 
important points concerning the general structure of the 
above flow equations. First, we note that t\a does not 
appear in the flow equations for the couplings /, g and 
u. This is because the flow equation for fi is given up to 
1-loop order by 
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Thus, the flow diagram is completely determined by the 
couplings /, g and u. This means that we do not need 



to know the expression of r\A- Another important point 
concerns the renormalization of 9. It is a known fact 
in topological field theory that the CS mass does not 
renormalize at all orders in perturbation theory and for 
the model we are considering it can be verified by explicit 
calculation up to 2-loop s Pq| . Thus, Zg = 1 and 9 = 
Za9q, which implies Eq.([l6[). 

The most important point concerning the above flow 
equations is related to the gauge dependence. It can be 
shown that the beta function for the gauge couplings are 
gauge independent in a minimal subtraction scheme [Q. 
However, is gauge dependent. When is evaluated 
at the fixed point it gives the r\ exponent. Note that 
exponents for the superconducting transition should be 
evaluated at the infrared stable charged fixed point, if wc 
assume that it exists. At the superconducting fixed point 
we must have tja = 1- This means that the fixed point 
value of a must be a = 0, the Landau gauge. Since at the 
neighborhood of the superconducting fixed point a flows 
to the Landau gauge, critical exponents are evaluated for 
a = and we shall fix it from the very beginning, as is 
customary in the literature. 

The explicit analytical expressions of i]^ and u are 
given by 
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The function u> and the anomalous dimension 77 have well 
defined limits as g — > 0, which correspond to the limit of 
the conventional GL model. As g — > 0, wc find r] — > — //4 
as g — > while lu — > 3n/2. The opposite limit, g — > 00, 
leads to the expected decoupling of the gauge and scalar 
degrees of freedom since both jn and u) tend to zero in 
this limit. The / 2 term in Eq.(|17j) comes from the graph 
shown in Fig.l. As mentioned in the introduction, this 
graph is zero when p — 0. In the scaling considered in 
this paper this graph survives since p ^ and we can 
recover the g — ► limit. 

The charged fixed points are given by /* = 16, g* 
arbitrary and 



l-2r)± W(2?7-l) 2 - 



160 



(23) 



where 77 = ??>/.(/*, 5*) and lu* = uj(g*). We have that u*± 
is real only for g* > g* fa 0.42, which corresponds to the 
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condition for the existence of charged fixed points in our 
model. The flow diagram is shown in Fig. 2 for g* = 0.48. 
The left charged fixed point, coresponding to u*_, is the 
so called tricritical fixed point. The tricritical fixed point 
is attractive over a line intercepting the origin called the 
tricritical line. The tricritical line separates the regions 
of first and second order phase transitions. The right 
charged fixed point, corresponding to u*,, is infrared at- 
tractive and describes the physics of second order phase 
transitions in superconductors. Note also in the flow di- 
agram the Gaussian and the XY fixed points. Another 
interesting point is that there exists another critical value 
°f 9*i 9c 2 ~ 0-81, such that for g* > g* 2 the tricritical 
fixed point is in the region of the plane uf defined by 
u < 0. In this region the tricritical fixed point lost its 
physical meaning. It results that the only charged fixed 
point is the infrared one. Consequently, for g* > g* 2 only 
second order behavior seems to be possible. 

Let us evaluate the critical exponents for the supercon- 
ducting transition in function of g* > g* . It is sufficient 
to evaluate rj and v since the other exponents are ob- 
tained from these ones via scaling relations. r\ is given 
simply by evaluating r]^ at the fixed point. The eval- 
uation of v, however, is more involved. In the critical 
theory this exponent is evaluated by considering a ip'tp 
insertion in the two point function. Thus, we must com- 

(1 2) 

pute the 1-particle irreducible function T\ ' subjected 
to the renormalization condition: 



rj' 2) (fei,fc 2 ; 



-k 1 -k 2 )\ 



S.P. 



1. 



(24) 



which will determine the renormalization constant of 
ip^ip, denoted by As usual the v exponent is 

given by the fixed point value of the RG function de- 
fined by 
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We find: 
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At 1-loop order we obtain 



The Fig. 3 shows a plot of the exponent v as a function 
of g* . Note that the plot is made for g* > g* ci which cor- 
responds to the region where charged fixed points should 
exist. We observe that as g* increases the value of v tends 
asymptotically to 0.6, which corresponds approximately 
to the XY value for the pure scalar model in the 1-loop 
approximation at fixed dimension. This is an expected 
result since for g — > oo the gauge modes decouple from 
the scalar modes. This result can be verified directly 
from the above RG functions by expanding for g large. 

In recent years it has been stablished that the v ex- 
ponent for the (non-topological) superconducting tran- 
sition is given by v = vxy ~ 0.67 |l{^||. Thus, we 
are tempted to search for what value of g* we have 
v = vxy- One obtains that v « 0.67 for g* 0.776. 
This is smaller than g* 2 and we have still the region of 
first order transition in the flow diagram. Unfortunately, 
a pathological behavior arises for this value of g* . The 
trouble comes from the rj exponent. Indeed, we have that 
rj w —2.47 < —1 while we know from the scaling relations 
that the condition 77 > —1 must be fulfilled. This result 
is probably an artifact of our 1-loop perturbation theory. 
Another related problem is that the above results fail in 
describing a flow diagram with charged fixed points in 
the limit g — > 0. In this limit our flow diagram resembles 
to that one obtained by Halperin et al. in the seventies 
Q and , as mentioned in the introduction, their result 
is an artifact of the e-expansion. Therefore, we should 
improve our perturbation theory in order to bypass all 
these difficulties. 



III. IMPROVED RG RESULTS 

As we have said in the last section, the beta functions 
for the topological GL model have a well defined g — ► 
limit. However, we have seen that in this limit no charged 
fixed points exist. Moreover, the 77 exponent could attain 
unphysical values which violate the scaling relations. The 
trouble is that there are in fact two fundamental lenght 
scales in this problem, namely, the correlation lenght £ 
and the magnetic field penetration depth, A. The £ is 
related to the scaling of the scalar field while A is related 
to the scaling of the gauge field. Thus, in principle, the 
renormalization conditions for the gauge coupling should 
be fixed at a different point of the scalar coupling. Since 
there is a relation between the lenghts £ and A, we must 
have also a relation between the corresponding renor- 
malization points. We develop this more general point of 
view by using a simple method suggested recently by Her- 
but and Tesanovic [111 . It consists in fixing the renormal- 
ization condition EqT(6) at the point given by k 2 — p 2 /c 2 , 
c giving in this way the ratio between the two scales of 
the problem. If we use the same reasoning here we obtain 
the following beta functions in the g — > limit: 
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where = — //4. Note that we have not exactly the 
same beta functions as in ref. jll| . This is due to the fact 
that we used the convention that all external momenta 
are incoming at the vertices. Our choice is more usual 
in the field theoretical literature and has the advantage 
that corresponding crossed graphs have the same value 
at the symmetrical point. However, these differences in 
the conventions does not change appreciably the physical 
results like the values of the critical exponents. 

The value of c can be fixed by demanding that, if 
charged fixed points do exist, it should happen at the 
same critical scale. At 1-loop this is equivalent to de- 
manding that the Ginzburg parameter n = u/2f should 
be invariant along the RG trajectoryconnecting the ori- 
gin and the tricritical fixed point pi[ . We choose the 
tricritical fixed point because there are good numerical 
estimates of k available |^,^l[ . If n tr i C denotes the value 
of k at the tricritical fixed point we obtain that c is de- 
termined by the following equation: 



c + 8- Vc 2 + 16c- 176 
40 ' 



(31) 



We take Ktric ~ 0.42/V2 which is the value obtained 
from Montecarlo calculations Thus, Eq.(|3l|) gives 

c « 27.78. For this value of c we have the flow diagram 
shown in Fig. 4 for the non-topological model. Fig. 5 
shows the detail of the region near the tricritical line 
which is not well seen in Fig. 4. 

The RG function in the g — > limit is given by 



1 




(32) 



We have the following result for the critical exponents: 



v w 0.676, 
r] w -0.14, 



(33) 
(34) 



a result in good agreement with the expected XY behav- 
ior. In ref. |ll| the value v w 0.53 was obtained by using 
a c corresponding to Ktric — 0.8/\/2 which is determined 
from the lattice dual model pi]| . However, they obtained 
v w 0.62 for Ktric = 0.42/V2- It is possible to improve 
this value obtained in ref. Jll[] by using directly Eq.(|25|) 
without expanding it, to obtain the better value v « 0.67, 
in close agreement with our result |22]j . This XY value 
has also be obtained in ref. It is worth to mention 
that non-perturbative RG calculations by Bergerhoff et 
al. ]l0| gives v w 0.53 in a certain approximation corre- 
sponding to a truncation of the average action (the Leg- 
endre transform of the Wilsonian effective action) at l^j 4 . 
Truncation at IV'I 8 gives the improved result v » 0.58. 



Let us come back now to the case of non- vanishing g. 
Since the g coupling is associated to the gauge degrees 
of freedom its momentum scale should be rescaled by c. 
Thus, in addition to Eqs.(|l7|) and (^) we have 



dq ( , cf 



(35) 



We keep the same value of c in order to preserve the 
features of the g — ► limit. Now we have charged fixed 
points even when g = and therefore there is no g* . On 
the other hand, we have still g* . We verified that g* 
is the same as before, that is, g* 2 w 0.81. This means 
that for g* > 0.81 we have only second order behavior. 
In Fig. 6 we show the flow diagram for g* = 0.48 with 
the detail of the tricritical region shown in Fig. 7. In 
Fig. 8 we plot the exponent v as a function of g* in 
the improved 1-loop calculation. Note that the shape 
of the plot is qualitatively the same as in Fig. 3. The 
main difference is that the vertical scale is compressed 
and the most remarkable feature is the fact that v is 
now much less sensible to the value of g* , showing only a 
small deviation from the XY scaling. For instance, the 
exponents v and r\ for g* — 0.48 are given by 



vk, 0.628, 
77 « -0.11. 



(36) 
(37) 



Of course, we have still the same g — > oo limit with expo- 
nents 77 = and v = 0.6, corresponding to the decoupled 
situation. 



IV. DISCUSSION 

The main aim of this paper is to initiate a careful study 
of a topological GL model from the point of view of criti- 
cal phenomena. For this reason we have concentrated the 
efforts on the RG flow and the evaluation of critical expo- 
nents. The results show that the topological coupling is 
a good control parameter with respect to the fixed point 
structure. For instance, we have seen that the region of 
first order transition is crunched as the topological cou- 
pling is increased and eventually the type I behavior is 
lost. An interesting point is that the v exponent does 
not fluctuate very much around the XY value. 

Besides superconductivity, the topological GL model 
may be useful in other physical contexts. For example, it 
can be applied in the study of soft materials like the chiral 
liquid crystals J23|. In this case the gauge field should 
be thought as a director field and the CS term is used 
in order to introduce the chirality for the constituting 
molecules. 

Another interesting problem is the physics of the chi- 
ral spin state. This state arises when one considers the 
Hamiltonian for the Heisenberg antiferromagnet whose 
spins interact not only through nearest neighbors in- 
teraction, but also through a next-nearest neighbor one 
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Q,^5[. The next-nearest neighbor interaction frustrate 
the Neel state and generates a mean field solution corre- 
sponding to a disordered spin state. The most stable con- 
figuration corresponds to the so called chiral spin state. 
The continuum effective theory is obtained by computing 
the fluctuations around such a mean field ground state 
and it contains a dynamically generated gauge field. The 
gauge sector of the effective action is identical to that 
one of the topological GL model. In this case the gauge 
couplings are given in terms of mean field parameters of 
the original model (for details see ref. p5[). 

Finally, we hope that this work will contribute to im- 
prove the understanding of GL models. 
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FIG. 1. Graph contributing to the f 2 term in the beta 
function for u. The dashed lines represent the gauge field 
propagator. 
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FIG. 2. Flow diagram for g* = 0.48. 



FIG. 5. Detail of Fig. 4 in the region near the tricritical line. 
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FIG. 3. Plot of v as a function of g* for g* > 0.42 



FIG. 6. Flow diagram for the topological model with im- 
proved 1-loop calculations for g* =0.48. 
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FIG. 4. Flow diagram for the non-topological model. 



FIG. 7. Detail of Fig. 6 in the region near the tricritical line. 
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FIG. 8. Plot of the exponent v as a function of g* . 
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